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The effect of nucleon-nucleon correlations in symmetric nuclear matter at finite temperature is 
studied beyond BCS theory. Starting from a Hartree-Fock description of nuclear matter with the 
Gogny effective interaction, we add correlations corresponding to the formation of preformed pairs 
and scattering states above the superfluid critical temperature within the in-medium T-matrix ap- 
proach, which is analogous to the Nozieres-Schmitt-Rink theory. We calculate the critical tempera- 
ture for a BEC superfluid of deuterons, of a BCS superfluid of nucleons, and in the crossover between 
these limits. The effect of the correlations on thermodynamic properties (equation of state, energy, 
entropy) and the liquid-gas phase transition is discussed. Our results show that nucleon-nucleon 
correlations beyond BCS play an important role for the properties of nuclear matter, especially in 
the low-density region. 

PACS numbers: 21.65.-f, 26.60.-c, 64.70.F- 



o 

=5 



> 

(N 

00 

in 

o 
o 



% 



I. INTRODUCTION 

Pairing and nucleon-nucleon correlations are impor- 
tant properties of interacting nuclear systems. For exam- 
ple, in the weak-coupling limit, i.e., at high density, the 
nucleons form Cooper pairs, and below a certain critical 
temperature T c the system is in a superfluid phase as de- 
scribed by the Bardeen-Cooper-Schricffcr (BCS) theory. 
In the strong-coupling limit, i.e., at low density, neutrons 
and protons form deuteron bound states which will con- 
dense if the temperature is below the critical temperature 
for the corresponding Bose-Einstein condensation (BEC). 
It was theoretically predicted .1] and recently confirmed 
by experiments with ultracold atomic Fermi gases @, Q 
that there is a smooth crossover between the BCS and 
BEC limits. Qualitatively, especially at zero tempera- 
ture, these features can be studied within the BCS (mean 
field) approximation [4J. Quantitatively, however, the 
critical temperature obtained in this way is too high be- 
cause the BCS theory does not include the existence of 
non-condensed pairs at finite temperature. In order to 
go beyond mean field, one has to consider pair correla- 
tions already above the critical temperature, as in the 
Nozieres-Schmitt-Rink (NSR) theory [lj. Especially in 
the low density region, where the coupling between nu- 
cleons is strong, such correlations modify the mean-field 
results to a large extent. 

At present, there are several groups who have stud- 
ied nuclear matter within the NSR approach. Pioneering 
work has been done by the Rostock group @, Q . There 
are also extensions where the correlations are considered 
in a more self-consistent way, like in the self-consistent 
Green's function method [7J, |8|. A generalization to tem- 
peratures below the superfluid transition temperature 
was discussed by Bozek |9j • In the case of ultracold Fermi 



gases, where the results can be compared with very pre- 
cise measurements, theories for the BEC-BCS crossover 
based on the NSR approach T0| have been very success- 
ful [HI. 

It is well known that there exists a liquid-gas phase 
transition in nuclear matter. Experimental information 
can be obtained from multifragmentation (see, e.g., [12| - 
Il5||). The critical temperature deduced from these ex- 
periments depends on the mass of the nuclei and can be 
as low as 6.7 MeV [15| in the case of small systems. For 
infinite nuclear matter, theoretical predictions give much 
higher values for the critical temperature between 14 and 
18 MeV [13, [H| (see Ref. QJ] for a recent theoretical 
study). Below that temperature, nuclear matter is unsta- 
ble in a certain range of low densities. Within mean-field 
theory, we know that the BCS-BEC crossover is com- 
pletely covered by the instability region of the liquid- 
gas phase transition. Nevertheless, the investigation of 
low-density nuclear matter is of interest for applications 
where regions of low density appear in the framework of 
the local-density approximation. Contrary to the nuclear 
matter case, the whole crossover can be studied in the 
case of ultracold atomic Fermi gases @, Q , because the 
pair correlations stabilize the gas [l| such that the system 
does not collapse into its solid ground state but it remains 
in its metastable gas state. By analogy, one expects that 
pair correlations will stabilize low-density nuclear matter 
and thus reduce the liquid-gas coexistence region. One 
of our subjects of investigation will be how strong this 
effect of nucleon-nucleon correlations on the liquid gas 
phase transition is quantitatively. 

Furthermore, in this paper we will calculate the equa- 
tion of state of hot and dense symmetric nuclear matter, 
taking into account the contribution of the mean field 
together with the nucleon-nucleon correlations. For the 



mean field we will use the Gogny interaction because it 
is known to give a good description of the single-particle 
and thermodynamic properties of nuclear matter, includ- 
ing saturation at the right density, the liquid-gas phase 
transition, etc. For the part beyond the mean field, we 
use the T-matrix (or ladder approximation) which con- 
tains the information on two-particle correlations. This 
also allows us to extract the critical temperature for pair 
condensation smoothly interpolating between the BEC 
and BCS regimes. 

The paper is organized as follows. In Sec. [Til we wm 
give a summary of the formalism. The numerical results 
are provided in Sec. IIIII The last section is devoted to 
the summary and discussions. 



II. FORMALISM 

Before explicitly including two-particle correlations, we 
calculate the single-particle Green's function within the 
Hartree-Fock (HF) approximation. In order to get a rea- 
sonable description of the single-particle energies, we use 
the density-dependent Dl Gogny effective interaction to 
describe the mean field. This force gives nuclear binding 
at the right saturation point and many other properties 
of nuclear matter and of finite nuclei [l7| ■ It has the form 



V(r) = J2 ( W ™ + B m P a - H m P T - M m P a P T )e- T ' 2 /^ 

m— 1 

+ t {l + x a P a )p a 6(r), (1) 

where the P a and P T are, respectively, the spin and 
isospin exchange operators. The spin-orbit coupling term 
is neglected here, since we consider only the properties 
of infinite nuclear matter. For the parameters we use the 
values given in Ref. [17j [34| . For details of the HF de- 
scription of nuclear matter at finite temperature with the 
Gogny force, see Refs. [HrHS]. The HF mean field Zhf 
contains the direct, the exchange, and the rearrangement 
contributions. Because of the finite range of the Gogny 
force, the exchange contribution is momentum depen- 
dent, and the single-particle Green's function takes the 
form 



G H f(p,u) 



where £ p is the quasiparticle energy defined by 



£ P = o ^hf(p) - li , 

Im 



(2) 



(3) 



where \i denotes the chemical potential. In order to facil- 
itate the numerical calculation of the correlation effects, 
we use the effective-mass approximation for the Gogny 
mean field, i.e., we write II 811 
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FIG. 1: The Feynman diagrams for the Green's function 
(top), for the self-energy (middle), and for the T-matrix in 
ladder approximation (bottom). 



There are different ways to define the effective nucleon 
mass m* . In principle, m* is momentum dependent 20]. 
Here we use the effective mass defined by expanding Eq. 
(|3]) around p = (we checked that for the final results it 
makes almost no difference if we expand around zero or 
around the Fermi momentum), i.e., 



J_ = J_ 2 dY> HF (p) 
m* m dp 2 



p=0 



fj, =/i 



^hf 



(0) 



(5) 
(6) 



However, the effective-mass approximation will only be 
used for the calculation of the correlation effects, while 
the mean-field contributions will be computed with the 
full momentum dependence of Yihf{p)- 

In principle, we are looking for the full single-particle 
Green's function G including correlations. The Dyson 
equation can be written as 



G- 1 (p,u>) = G^ l F (p,u>)-£( P ,u>) 



(7) 



where S is the correlation contribution to the single- 
particle self-energy. Since the Gogny force is a density- 
dependent effective interaction, which is designed to give 
good results already at the HF level, we suppose that 
the Gogny mean field accounts already for most of the 
correlation effects. We therefore demand that the cor- 
relations do not shift the quasiparticle energies £ p , i.e., 
^(Pi£p) — 0j ancl that the role of the correlations is just 
to reduce the strength of the quasiparticle pole and to 
distribute the remaining strength in the continuum of 
the spectral function. Hence, we define £ to be the self- 
energy subtracted at £ p : 



E(p,w) = E(p,w) -Re£(p,£ p 



(8) 



In order to describe pair correlations, we calculate the 
self-energy S within the T-matrix or ladder approxima- 
tion, as shown in the lower part of Fig. [H This is a fre- 
quently used lowest-order correction [l], 0, [(| [I(| , imply- 
ing, however, that vertex corrections as well as screening 
of the interaction due to the medium effects are neglected. 
Since our aim is not a completely self-consistent de- 
scription of the spectral function as in the self-consistent 



Green's function method 0, [8j , we make the assumption 
that the correlations can be treated as a small correction 
to the Gogny HF self-energy. This allows us to use the 
HF Green's function Ghf in the calculation of the T ma- 
trix and of the self-energy E. Then, for consistency, one 
should also keep only the first-order term of Eq. (O, i.e., 



G(p, u) = G HF (p, u) + G 2 HF (p, w)S(p, w) 



(9) 



A diagrammatic representation is given in the upper part 
of Fig. HI 

That the self-energy in T-matrix approximation should 
only be treated in first-order perturbation theory may 
also have a more formal reason. The T-matrix approx- 
imation corresponds to particle-particle random-phase 
approximation (pp-RPA) [21(. It can be shown that the 
ground-state energy calculated from the single-particle 
Green's function with self-energy in first order and in T- 
matrix approximation yields exactly the pp-RPA ground- 
state energy (22|. At least this holds true for the self- 
energy without subtraction procedure. Therefore our 
formalism is closely related to that of Ref. |23J], where 
the pp-RPA formalism is used, except that we apply the 



subtraction prescription while the authors of Ref. 23 1 
are obliged to reduce the correlation contribution by in- 
troducing a cutoff and to change the parameters of the 
Gogny force in order maintain the right saturation point 
of nuclear matter. 

Note that our approximations are analogous to NSR 
theory [l[, except that in NSR theory free Green's func- 
tions instead of HF ones are used and consequently no 
subtraction is made in the self-energy. In the case of 
nuclear matter, however, we cannot expect to obtain a 
good description of the full self-energy from such a simple 
model for the T matrix. This is why we use the Gogny 
mean field and the subtraction method described above, 
while the subtracted self-energy serves only to provide 
the energy dependence corresponding to the pair corre- 
lations in the channels we want to study. 

In order to get a simple expression for the T matrix, 
we use the separable Yamaguchi potential [24| , 



V a (k,k') = -X a v(k)v(k') 



(10) 



where k and k' are the incoming and outgoing relative 
momenta in the center-of-mass frame, and the form fac- 
tor is given by 



v(k) 



1 



k 2 +(3 2 



(11) 



As in Ref. |6|, we consider only S-wave scattering (a = 
1 Sq, 3 Si) and neglect the coupling between the 3 Si and 
3 -Di channels (which comes from the tensor force). With 
the parameters f3 = 1.4488 far 1 , Ai So = 2994 MeV far 1 
and Xs Sl = 4264 MeV far 1 Q, the low-energy nucleon- 
nucleon phase shifts and the vacuum binding energy of 
the deuteron (E® = —2.225 MeV) are very well repro- 
duced, see results for n = in Figs. [5] and [3J so that it 
is unlikely that the coupling between the 3 £>i and 3 Si 
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FIG. 2: The deuteron binding energy in nuclear matter from 
the Yamaguchi potential and including the effect of the Gogny 
mean field, as a function of the density for different values of 
the deuteron momentum K. The temperature is T — 10 MeV. 



channels would strongly modify our results. With the 
separable interaction, the resummation of the ladder di- 
agrams shown in the lower part of Fig. [T] reduces to a 
simple geometrical series, and the T matrix can be writ- 
ten as 



T a {k,k' ) K > u) = 



Va(k,k') 

l-J a {K,uj) ' 



(12) 



where k and k' are the incoming and outgoing momenta 
in the center of mass frame, K is the total momentum, 
and 



J a (K,u) 



d 3 k 
(2^)3 



V a (k,k) 



1 — /(^K/2+k) ~ /(£k/2-1c) ,,„x 

w _ ^K/2+k - ^K/2-k + «0 



The function /(£) = l/(e«/ T + l) is the Fermi function, T 
being the temperature. Within the effective mass approx- 
imation, Eq. (0|, the denominator of Eq. (fTB"]) does not 
depend on the angle between k and K, and the angular 
integral can be done analytically. The main contribution 
to the integral over the relative momentum comes from 
low momenta due to the form factor of the Yamaguchi 
interaction (k < 0). 

In the 3 Si channel, it can happen that J3s x {K, Wb) = 1 
at some energy uji, below the threshold energy 



to (K) 



K 2 

4m* 



2/x* 



(14) 



This means that there is a bound state (the deuteron) 
with binding energy Et,(K) — u>b(K) — ujq(K). As an 
example, the deuteron binding energies for different val- 
ues of the deuteron momentum K are displayed in Fig. 
[5] As one can see, the binding gets weaker with increas- 
ing density, and eventually the deuteron gets unbound 
at the so-called Mott density. Since the Pauli blocking 
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FIG. 3: In-medium scattering phase shift in the Si channel 
for K = as a function of E = fc 2 /m* for different densities 
and T = 10 MeV. 



effect gets weaker with higher deuteron momentum K, 
there exists for any density a Mott momentum Kmou 
above which the deuteron stays bound. 

The in-medium nucleon-nucleon phase shifts 5 a can 
easily be obtained from 1/(1 — J a ) = e t<5 °/|l — J a \. As 
an example, we show in Fig. [3] the phase shift in the 
3 Si channel for K — at different densities, as function 
of the energy E = lu + 2/j,* = k 2 /m*. We see that at 
higher densities, e.g., at n > no/5 (no — 0.17 fm~ 3 being 
the saturation density of nuclear matter), the phase shift 
is negative in the low-energy region and then becomes 
positive as the energy increases. The energy where the 
phase shift crosses zero is w = 0, i.e., E = 2/i*. At 
lower densities, when jj,* is negative, the phase shift is 
positive at low energy. At some very low density, the 
value of the phase shift at E — changes from to n. 
This happens precisely at the density below which the 
deuteron is bound. 

In terms of the T matrix, we can write the self-energy 
E depicted in the middle of Fig. Q] within the Matsubara 
formalism as 



E(p,iw„) = | 



E *T 



d 3 P ' 

(2tt) 3 



GHF{p',iUn>) 



xT a (k,k,K,iu n + iu n >), (15) 



where w„ and u n i are Fermionic Matsubara frequencies 
[u„ = (2n + l)7rT],k= (p-p')/2, andK = p + p'. The 
factor 3/2 is the product of a factor 1/4 from the aver- 
aging over spin and isospin in symmetric nuclear matter, 
of a factor 2 from the sum of direct and exchange contri- 
butions, and of a factor (25+ 1)(2T+ 1) = 3 for a = 3 Si 
and 1 So from the sum over spin and isospin in the loop. 
Using standard techniques [25j , the self-energy can be an- 
alytically continued to real energies, which is necessary 
for the calculation of the subtraction term E(p, £ p ) in Eq. 
©■ 

Inserting the self-energy into Eq. ^ , we calculate the 



density from 

n(T, a) 



-4rW 



_fk_ 



G(k,iu!r, 



(16) 



The factor 4 comes from the sum over spin and isospin. 
It is clear that the first term of Eq. © just gives the 
Hartree-Fock density, and the second term gives the cor- 
rection beyond the mean field approximation. After a 
lengthy derivation (see Appendix), one finds the follow- 
ing formulas initially given in Refs. [5j, [6j: 



n HF 



n HF + ribo 



The bound-state contribution reads 

d 3 K 



Inbound O 



g{ub{K)) , 



(17) 



(18) 



ik>k mm (2?r) 3 

where g(uj) = l/(e UJ ^ T — 1) is the Bose function. This 
term gives the nucleon density corresponding to a Bose 
gas of deuterons. The scattering-state contribution reads 



ll'scatt 



-6 



K>K A 



-6 



E 

3 Si, 



So 



d 3 K 
(2^)3 
d 3 K 

(2tt) 3 



duo 
u {K) * 

x (S a - 



g(u) 



did 

sin 25, 



(19) 



In Ref. [5j, these equations were derived in a different 
way using the optical theorem, analogously to the deriva- 
tion of a similar formula for the electron-hole system in 
Ref. [26|. 

Note that in spite of the double pole of the derivative 
of the Bose function at u) = 0, the integrand in Eq. (jT^j) 
has no pole. This is because 8 a crosses zero at u = 0. 
This simple zero is raised to a double one due to the 
difference of the two terms in the second line of Eq. (|T9]) 

Once we have calculated the density, we can calculate 
the pressure. To that end, we integrate the thermody- 
namic relation n — (dP/ da)^ over /x, i.e., 



p ( T > M) = / n ( T > v')d\x . 



(20) 



Then we calculate the free-energy density F/V, the en- 
tropy density S/V, and the energy density E/V from the 
thermodynamic relations 



F = -PV + fmV, S=- 



dF 
df 



and E = F+TS. 
(21) 



III. NUMERICAL RESULTS 

A. Density and the super-fluid critical temperature 

We calculate the total density by numerically integrat- 
ing Eqs. (HHJ) and (|T9")) . The results for the densities 
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FIG. 4: The densities at T = 20, 15.9, 10, and 5 MeV 
(from left to right) as functions of the chemical potential 
within Gogny HF (dashes) and with correlations (solid line). 
T liq-gas _ 15 g Me y j g the critical temperature for the liquid- 
gas phase transition. 



at different temperatures as functions of the chemical 
potential 



are shown in Fig. 21 Comparing the re- 
sults with correlations (solid lines) with the Gogny HF 
results (dashed lines), one can see that, for a given chem- 
ical potential, the correlations increase the densities. In 
the high-density region, we notice that the results with 
and without correlations converge to the same value, i.e., 
the correlations fade away at high density, as this can be 
expected. For example, at T = 5 MeV, the two results 
coincide starting from n — 0.07 fm~ 3 . This is a conse- 
quence of the Mott mechanism, which has been discussed 
at length in Ref. [5|. As mentioned above, the critical 
number density where the bound state (at K — 0) disap- 
pears is called Mott density. When we change the tem- 
perature from 5 MeV to 10, 15.9, and 20 MeV, the Mott 
density changes from 0.07 far 3 to 0.12, 0.18, and 0.22 
fm~ 3 . This means that the mean field approximation is 
valid in the high density region. Below this region, the 
contribution of the nucleon-nucleon correlations is impor- 
tant. 

From this figure we also can see that when the tem- 
perature is less than some critical value (T^ q ~ gas = 15.9 
MeV), the number density has three values correspond- 
ing to one definite value of chemical potential. This is a 
typical feature of the liquid-gas phase transition in nu- 
clear matter. We will discuss this phenomenon in detail 
in the next subsection. 

To see how large the correlation contribution to the 
density is, we show the composition of the system at dif- 
ferent temperatures in Figs. [5] and [51 Since the density 
ratios are shown as functions of the density and not of 
the chemical potential, there are unique solutions even 
for temperatures below Tl %q ~ 9as . In Fig. [5] one can see 
that at T = 5 MeV the correlation contribution to the 
total density is important at low density (n < no/4). At 
n = 0.02 fm~ 3 , the correlated part is even larger than 
the HF part. This means that most of the nucleons are 
in correlated pairs in this density region. With increas- 
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ing temperature, e.g., at T — 10 MeV as shown in Fig. 
|6j the ratio of the correlated density to the total den- 
sity decreases, but the density region with sizeable nu- 
cleon correlations is enlarged. Here we do not separate 
the correlation contribution into bound and scattering 
state contributions, since individually they are not very 
meaningful, as discussed in Ref. [5J. For instance, if the 
temperature is much higher than the deuteron binding 
energy, the first term of the scattering-state contribution 
(IT91 cancels almost exactly the bound-state contribution 

dm. 

In the above calculation, when the temperature is be- 
low some critical value, we get a divergence in the T ma- 
trix. This pole corresponds to the formation of Cooper 
pairs at high density and to Bose-Einstein condensation 
of deuterons at low density. Below this critical tempera- 
ture T c , the equations for the density of the system are 
not applicable any more. In the superfluid phase, one 
would have to include the nucleon pairing gap explic- 
itly in the single-particle Green's function (which then 
becomes a 2 x 2 matrix in Nambu-Gorkov space [25j). 
which is beyond the scope of this paper. However, we 
can determine the critical temperature of the superfluid 
transition as the temperature where the T matrix devel- 




n (fm ) 



FIG. 7: Superfluid critical temperature as a function of the 
(total) density. The solid line is the full calculation, while the 
long dashes correspond to the BCS result. The short dashes 
show the critical temperature of Bose-Einstein condensation 
of a deuteron gas. 



ops a pole at zero total momentum (K = 0) and at zero 
energy (u> = 0). This is the well- known Thouless crite- 
rion [27j | for the onset of superfluidity, coinciding with 
the BCS gap equation when the gap A goes to zero: 



1- J a (K = 0,u = 0;T = T c ) 



0. 



(22) 



From this equation we get the critical temperature as 
a function of the effective chemical potential. Using 
the relation between the effective chemical potential and 
the number density, we obtain the superfluid region be- 
yond the BCS (mean field) result as shown in Fig. [7] 
Qualitatively, this result is similar to the one in j6| ex- 
cept that we have a lower critical temperature for the 
superfluid phase transition. The maximum T c in [6| is 
7.2 MeV at n = 0.12 fur 3 , while we have T c = 4.5 MeV 
at n = 0.05 fm -3 . The difference stems from the Gogny 
mean field, in particular from the effective mass, which 
was neglected in Ref. [g. One realizes that a T c of 4.5 
MeV is still very high, leading to a maximal gap of about 
7 MeV, about three times as much as the maximum value 
of the neutron-neutron gap in the spin singlet channel. 
The reason clearly stems from the slightly stronger at- 
traction in the proton-neutron isoscalar channel. How- 
ever, in finite nuclei barely any enhancement of pairing 
in the S = 1, T = channel can be detected. Proba- 
bly important screening is at work in that channel. In 
nuclear matter, this has been investigated in Ref. [31j . 
The addition of screening effects is, however, beyond the 
scope of this paper. 

As mentioned above, T c as a function of fi coincides 
with the BCS result. As a function of the density, the 
difference between the results T c (n) with and without 
correlations comes only from the different relations for 
n as a function of /i. Since the correlation contribu- 
tion to the density vanishes in the high density region, 
the phase boundary coincides with the BCS curve (long 
dashed line, which is obtained with uhf only). At very 
low density and temperature, the main contribution to 



the density comes from the deuteron bound state, as can 
be seen from Eqs. (fTT]) and ([18]), Close to the Bose 
critical temperature, the Bose distribution function in 
Eq. (|18j)) starts to diverge and, therefore, dominates the 
whole expression for the density. Therefore the super- 
fluid critical temperature at low density coincides with 
the critical temperature for Bose-Einstein condensation 
of a deuteron gas, which is given by 



T. 



BEC 



m \6C(3/2) 



2/3 



(23) 



(with £(3/2) = 2.612 . . . ) and is shown as the short- 
dashed line in Fig. [7] 

A surprising behavior of our result is that in the density 
region between 0.04 fm~ 3 and 0.05 fm~ 3 , Eq. ([22]) for 
the critical temperature has three solutions for one given 
density. This behavior is not easy to understand from 
physical intuition. It seems to be_related to the effective 
mass, since it is absent in Ref. [6]. Anyway, as we will 
show in the next subsection, this density region lies inside 
the unstable region of the liquid-gas phase transition. 



B. Pressure and liquid-gas transition 

As it was shown in Fig. [4] there is a region of den- 
sities where the chemical potential decreases with in- 
creasing density. This is a typical feature of a liquid-gas 
phase transition. In order to determine the boundary of 
this first-order phase transition, we need the pressure. 
In principle, one can get the pressure as a function of 
temperature and chemical potential, P(T,(j,), from the 
number density n(T, fj,) by integration over the chemical 
potential /i, cf. Eq. (|20[) . However, since there is a first- 
order phase transition, n is not a single- valued function 
of [i any more. We therefore transform the integral over 
/i into an integral over Uhf' 



P(T,n HF ) = 



n HF {T,fi) 



n(T,n' HF ) 



dji 



dn 



HF 



dn'. 



HF- 



(24) 



Since \i is a single-valued function of uhf (see dashed 
line in Fig. 0]), this integral is well defined. In this way 
we obtain the pressure as a function of uhf, but neither 
Uhf nor P are single-valued functions of \x. 

If we plot the pressure as a function of the total 
density n instead of uhf, we get the results shown in 
Fig. [8] Unfortunately, we cannot calculate the pressure 
for T < 4.5 MeV, at least not at densities above 0.05 
fm~ 3 , because our method to calculate the pressure at a 
given density n necessitates the calculation of all densi- 
ties n' < n, i.e., including the density at n = 0.05 fm -3 
where T c is maximum. For comparison, we also give the 
results for the pressure within the mean-field approxima- 
tion (dashed lines in Fig. [8]). As it can be seen, the 
main effect of the nucleon-nucleon correlations is to in- 
crease the pressure at very low densities. However, in the 
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FIG. 8: The pressure as a function of density at different 
temperatures. Solid lines: with correlations; dashed lines: 
mean-field results. 



16 


1 

- without correlation y-'yr^s 

with correlation S yf \v\ 


- 


14 


/ / 1 »\ \ 


- 


12 
10 


'' /' / \ \ 
/ / / \ \ 

— ■** / ' / * 


- 


8 
6 


/ / t 1 

^^ >' / 

/ 1 / 


i- 



0.001 



0.01 
n (fm~ 3 ) 



0.1 



FIG. 9: The liquid-gas phase diagram as function of density 
and temperature (for T > 4.5 MeV). The thin lines are the 
boundary of the coexistence region, while the thick lines are 
the boundary of the spinodal region. Solid lines: with corre- 
lations; dashed lines: mean field results. 



case T = 5 MeV shown in Fig. [H the pressure at higher 
densities is lower than the HF result. 

Using the pressure, one can determine the coexistence 
region of the liquid and gas phases of nuclear matter from 
the following conditions: 

P(T,m) = P(T,n 2 ) and n{T, ni ) = m(2> 2 ). (25) 

The result is shown in Fig. [5] as the thin solid line. At 
the same time, we can determine the spinodal curve from 
the zeros of dP/dn (or, equivalently, of 3/i/<9n), which 
is shown as the thick solid line in Fig. |9] In the re- 
gion under the spinodal curve, the system cannot exist 
in a homogeneous phase. In the region between the thin 
solid line and the spinodal curve, the gas phase (left-hand 
part) or the liquid phase (right-hand part) can exist as 
a metastable state. For comparison, the corresponding 
mean-field results arc presented in Fig. [9] as the dashed 
lines, which coincide with Fig. 6 of Ref. [20]. 

Comparing the results with and without correlations, 
one can see that the correlations decrease the phase- 
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FIG. 10: Phase diagram combining the boundary of the su- 
perfluid phase (lower curve), the liquid-gas coexistence region 
(upper curve), and the spinodal line (dashed curve). The rea- 
son why the spinodal and coexistence curves end at T c and 
4.5 MeV, respectively, is not physical but it simply means 
that our model does not allow us to compute them at lower 
temperatures (see text). 



transition temperature in the low-density region and re- 
duce the unstable region of the liquid-gas phase transition 
considerably. As mentioned in the introduction, this is 
an expected result. In the high density region, the effect 
of the correlations is almost negligible. 

We can determine the critical temperature of the 
liquid-gas transition, i.e., the maximum temperature of 
the coexistence and the spinodal curves, from 



dP 
dn 



8 2 P 



.'jitt-tj 



dn 2 



= 0, 



(26) 



see Fig. M In this way, we obtain T^s as = 15.9 MeV, 
which coincides with the mean- field result [19|, [2CJ . The 
fact that T^ q ~ 9as remains unchanged is an artifact of our 
present approach to treat the correlation effects only at 
a pcrturbative level, as explained in Sec. HU As shown in 
Ref. [30], the inclusion of dcutcron (and heavier) clusters 
should reduce the liquid-gas critical temperature. We 
would have to do the calculation more self-consistently in 
order to get a lower critical temperature than the mean- 
field result. 

In Fig. [lOl the results of Fig. [7] for the superfluid crit- 
ical temperature T c (lower solid line) and Fig. [9] for the 
liquid-gas coexistence region (upper solid line) and the 
spinodal instability region (dashed line) have been com- 
bined in a single phase diagram. As explained above, we 
unfortunately cannot calculate the liquid-gas coexistence 
curve for T < 4.5 MeV, but extrapolating the solid curve 
to lower temperatures and remembering that at T = 
the liquid phase gets stable at saturation density, it is 
clear that the coexistence curve will cross the superfluid 
T c curve at n ~ no, i.e., as one would expect, homoge- 
neous nuclear matter with pairing is stable above this 
density. From the results of Ref. [28| one can presume 
that the liquid-gas coexistence region will be slightly re- 
duced below the superfluid critical temperature, but this 
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FIG. 11: Energy per nucleon as a function of density for dif- 
ferent temperatures. Solid lines: with correlations; dashed 
lines: mean field results. 
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FIG. 12: Entropy per nucleon as a function of density for 
different temperatures. Solid lines: with correlations; dashed 
lines: mean field results. 



effect should be almost negligible in the case of symmetric 
nuclear matter considered here |28|,[29|. At low densities, 
superfluid matter is never stable, because the superfluid 
T c curve stays always below the coexistence curve. 

The spinodal curve (dashed line) can be calculated un- 
til it reaches the superfluid region. From this we see that 
superfluid nuclear matter is metastable below n ~ 0.045 
fm~ 3 and above n ~ 0.1 fm~ 3 . Note that on the low- 
density side, the density region where the gas phase is 
metastable is strongly increased by the correlations, es- 
pecially when we approach the superfluid transition tem- 
perature. This confirms our expectation mentioned in 
the introduction that the correlations have a stabilizing 
effect. However, the BEC-BCS crossover lies still in the 
unstable region of the liquid-gas phase transition. 



C. Energy and entropy 

The energy and the entropy can be obtained from the 
pressure with the help of the thermodynamic relations 
(|21|). Results for the energy per nucleon, E/A, and for 
the entropy per nucleon, S/A, for different temperatures 
are shown in Figs. 1111 and [T2l The corresponding mean- 
field results (dashed lines) are also shown for compar- 
ison. The results shown in Fig. [TT] indicate that, for 
fixed temperature, the correlations shift the minimum 
of the energy per nucleon to slightly higher densities. 
Fortunately the change is very small, because otherwise 
we would have to readjust the parameters of the Gogny 
force, which gives the right saturation density and energy 
at zero temperature without correlations. 

In the low-density region, where the dcutcrons and the 
nucleon- nucleon scattering states dominate, the energy 
per nucleon is lower than that the HF result. When the 
density is high, the correlation effect goes to zero and the 
energy per nucleon gets close to the mean-field result. 

When the density approaches zero, both results go to 



the classical value of an ideal gas of nucleons, 

lim E/A = §T. (27) 

This is not surprising, since even the lowest temperature 
considered here, T — 5 MeV, is still much higher than 
the deuteron binding energy so that almost all deuterons 
will be dissociated. However, the result ([27]) is also found 
at temperatures much lower than the deuteron bind- 
ing energy. This is because, at finite temperature, the 
deuterons are always dissolved in the low-density limit. 
This is a consequence of the mass-action law and can eas- 
ily be understood as follows: At low density, the chem- 
ical potential of the nucleons, /i, gets strongly negative, 
fi <C — T. The chemical potential of the deuterons is 
2/j,, which is even more negative. So the nucleon density 



n/T 



is much larger than the deuteron density oc e 



2/i/T 



oc e 

Only at zero temperature, where the system remains a 
deuteron BEC at arbitrarily low densities, the energy per 
nucleon approaches —1.12 McV (half the deuteron bind- 
ing energy) in the limit n — >■ Q. 

The results for the entropy (cf. Fig. fTS]) have been 
calculated from Eq. (|2"Tj) and show that, for fixed tem- 
perature, the entropy per nucleon decreases with increas- 
ing density. In the zero-density limit, the entropy per 
nucleon increases logarithmically, in agreement with the 
result for a classical ideal nucleon gas. As is clear from 
the discussion above, the correlations do not change this 
asymptotic behavior. At slightly larger values of the den- 
sity, the correlations tend to increase the entropy. 



IV. SUMMARY 

In this paper, we discussed the effect of pair correla- 
tions beyond the mean-field approximation in symmetric 
nuclear matter above the superfluid critical temperature. 
We include the effects of non-condensed pairs (deuterons) 
as well as the contribution of scattering states. For the 



mean field, we use the Gogny effective interaction in order 
to get the right saturation properties of nuclear matter. 

Starting from the single-particle Green's function 
within the Gogny HF approximation, we include the cor- 
relations in a perturbative way by considering in addition 
to the HF Green's function the diagram with one self- 
energy insertion, the self-energy being calculated in lad- 
der approximation. This approximation scheme is analo- 
gous to the well-known NSR approach. However, in order 
to avoid double counting of the quasiparticle energy shift 
which is already accounted for by the Gogny mean field, 
we have to subtract the self-energy at the quasiparticle 
energy. This leads finally to the same formula for the 
density in terms of the in-medium scattering phase shifts 
as given in Ref. [a]- We use a separable Yamaguchi 
potential in order to get an analytical formula for the 
in-medium T matrix and the phase shifts. 

Evaluating numerically these formulas for the den- 
sity, we discussed the different density contributions 
in hot and dense nuclear matter and found that the 
nucleon-nucleon correlations are important in the low- 
temperature and low-density region (n < hq). The cor- 
relation effect on the superfluid critical temperature was 
discussed. The result interpolates between the critical 
temperature for Bose-Einstein condensation at low den- 
sity and the BCS critical temperature at high density. 
We found that the maximum of the superfluid critical 
temperature decreases from 7.2 MeV (the value given in 
Ref. @) to 4.5 MeV when the effective mass to* due to 
the Gogny mean field is taken into account. 

Then we studied the liquid-gas phase transition in hot 
and dense nuclear matter with the help of the pressure 
calculated from the density. Especially at low density, we 
found that the boundaries of the coexistence and spinodal 
regions of the phase transition are shifted by the pair cor- 
relations. As we expected, the stable and metastable re- 
gions of the gas phase are strongly enlarged. In particular 
near the superfluid transition temperature, the gas phase 
stays metastable up to much higher densities if the corre- 
lations are taken into account. However, the correlations 
are not strong enough to suppress the liquid-gas transi- 
tion. This could have been anticipated from the fact that 
the liquid-gas critical temperature is much higher than 
the superfluid one [28|, [29j . Because of our perturbative 
treatment of the correlations, the critical temperature of 
the liquid-gas transition remains the same as within the 
mean field approximation. 

Finally, we calculated the energy and entropy of 
nuclear matter from thermodynamic relations. The 
nucleon-nucleon correlations decrease the energy per nu- 
cleon in the low density region but increase it at high 
density. For the entropy, the correlations always give a 
positive contribution. 

As mentioned before, our result for the critical temper- 
ature of liquid-gas phase transition is not affected by the 
pair correlations because they are treated only perturba- 
tively. One should improve this by taking the correlations 
into account self-consistently. Then the correlations will 



have some effect on the HF field and the critical tem- 
perature will change. The saturation point of nuclear 
matter, given correctly by the Gogny interaction within 
the HF approximation, may be changed, necessitating a 
readjustment of the parameters of the Gogny force. 

Our equation of state is only valid for temperatures 
and densities above the superfluid critical temperature. 
In order to get a result which is valid in the whole temper- 
ature and density plane, one should introduce the pairing 
gap A into the single-particle Green's functions. Some 
work in this direction has been done for nuclear matter 
9] , and quite elaborate theories have been developed for 
the BEC-BCS crossover in ultracold atomic Fermi gases 
|10( . We leave this for future study. Another important 
extension of the present work would be to consider the 
case of asymmetric nuclear matter and neutron matter, 
since these are of great importance for the study of neu- 
tron stars and their formation. 
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Appendix: Derivation of the density formulas 

In this Appendix we give a more transparent derivation 
of the density formulas (fTTl) . (TT51) . and p^|) . which were 
initially derived in Refs. 5, 26]. For better readability, 
we will not write out the sum over a — Si, 1 Sq and 
suppress the index a in this appendix. 

Let us recall the spectral representation of the T ma- 
trix, 



T(k,k',K,uj) = V(k,k')- 



du'ImT{k,k',K,uj') 

IT LO — bj' + i0 



where uj can be real or complex. Analogous dispersion 
relations exist for the self-energy £(p, w) and for the two- 
particle propagator J(K,tu) defined in Eq. (fl3|) . Using 
Eq. (|28[) . one can evaluate the frequency sum in Eq. 
(|15p. and one obtains the well-known expression for the 
imaginary part of the self-energy: 

ImS(p,w) = - / — — ^ImT(k,k,K,u + £ p >) 

x [/&,') + <?(w +&>')]• (29) 

where k and K are the relative and total momenta as 
defined below Eq. (fT5j) . 
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The correlation correction to the density is given by 
n corr = -4T2^ / 77^3 — — ~ ■ ( 30 ) 



(2tt) 3 



(w„ - £p) 2 



If we use the spectral representation of £, the frequency 
sum can be evaluated with the result 



= -4 



d 3 P vf^i w ^IM_IM 



(w ~ £p) 5 



(31) 

where T 3 denotes the principal value. Inserting Eq. (|2"5)l 
into this expression, one obtains with the help of the 
relation /(&)/(&,) = <?(£ P + &0[1 - /(&) - /&/)] and 
after some transformations 



d 3 p d 3 p' 

(2tt) 6 



— ImT(k,k,K,u}) 



x(i-/fc)-^) tf7 ( vy - ^ 

l w Sp sp' J 

The next step is to introduce the new variable u/ = £ p + 
£ P ' and to replace the integral over p' by an integral over 
ijj' . Then, using the imaginary parts of Eqs. (|12p and 
(|13p . one can show that the resulting expression for n corr 
can be rewritten as 



= 6 



d 3 K fdujduj' 
VI — Im 



(2tt 



1 



1- J(if,w) 



xlmJ^cO ^-y . (33) 
With the help of the dispersion relations for the real 



parts, this expression can be further reduced to 
f d 3 K f du . . / 1 d 

ncorr = 6 J (2^F J ~ g{uj) { T—j^ ReJ 



T T d 1 

Im J — — Re ■ 



did 1 — ,/ 



(34) 



(the arguments of J(K, u>) have been suppressed for 
brevity) . In order to express everything in terms of the 
in-medium scattering phase shifts S = — Imln(l — J), we 
notice that 

dS 1 d ^ ^ 1 d , _, 

— =Im- -— ReJ + Re -— ImJ, (35) 

aw 1 — J duj 1 — J au 



= sin 5 cos 5 — | sin 2(5 . (36) 



Im J Re 



With these relations, Eq. (|34|) can be rewritten as 

„ f d 3 K f duj , d , , . , 
" corr=6 i (2^7 V ff(w) AJ^ 2 sm2(5 ) • (37) 

The final step is to integrate by parts over u and to 
separate in the resulting integral the contributions of 
uj > ujq(K) (scattering-state contribution n SC att) and 
uj < loq(K) (bound-state contribution ribound)- The latter 
reduces to Eq. ((18)) since the phase shift below threshold 
is (see also Fig. 7 of Ref. [l|) 



6(K,u><uo(K)) = 



if K < K 



Motti 



o, 

nO(u-u b (K)), iiK>K Mo tt- 

(38) 
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